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TECHNICAL  NOTE  I9O3 

EXPERIMENTAL  INVESTIGATION  OF  MOVING  PRESSURE 
DISTURBANCES  AND  SHOCK  WAVES  AND  CORRELATION 
WITH  0NE-4)IMENSI0NAL  UNSTEADY-FLOW  THEORY 
By  Paul  W.  Huber,  Cliff  E.  Fitton,  Jr.,  and  F.  Delpino 


SUMMARY 


Experimental  pressure-time  measurements,  along  with  sohlieren 
photographs  of  the  unsteady— flow  phenomena  initiated  in  constant— cross- 
sectional-area  passages,  are  presented.  The  methods  of  calculation  of 
the  flow  phenomena  are  based  on  the  assumptions  of  one-dimensional, 
nonviscous,  unsteady  flow.  The  unsteady  flow  —  which  is  studied  both 
experimentally  and  theoretically  —  includes  propagation  of  pressure 
disturbances,  which  may  be  either  expansion  zones  or  discontinous 
compression  fronts,  temperature  contact  discontinuities  in  the  flow, 
and  the  reflections  and  interactions  of  these  flow  quantities.  A  very 
close  agreement  is  shown  between  the  experimental  measurements  of  the 
flow  and  the  flow  predicted  by  the  theory  over  a  wide  range  of  flow 
conditions.  Special  instrumentation  was  developed  during  these  tests 
which  greatly  facilitated  the  measurement  of  the  unsteady-flow  phenomena. 


* 

INTRODUCTION 


The  important  role  played  by  moving  pressure  disturbances  in 
various  aerodynamic  configurations  is  becoming  increasingly  apparent. 

In  pulse— jet  engines  and  heat  engines,  such  as  the  comprex  or  pressure 
exchanger  (reference  l),  moving  compression  and  expansion  waves  are 
alternately  used  to  compress  the  gas  before  heat  is  ad.ded  and  to  expand 
the  hot  gases  in  the  same  passage,  in  much  the  same  manner  as  is  obtained 
by  the  action  of  a  mechanical  piston.  The  action  of  pressure  disturbances 
in  the  subsonic  part  of  a  supersonic  diffuser  is  recognized  as  causing  an 
instability  of  the  shock  wave  with  a  resulting  adverse  effect  on  pressure- 
recovery  efficiency.  In  some  crises  (reference  2),  constant— geometry 
supersonic  diffusers  have  been  started  at  Mach  numbers  below  the  theoreti¬ 
cal  starting  Mach  number  (based  on  steady— flow  considerations)  by  means 
of  pressure  disturbances.  In  steady— flow— type  burners,  as  well  as  in 
combustion  engines,  the  action  of  pressure  disturbances  and  detonation 
waves  produce  a  marked  change  in  the  combustion  phenomena. 


Preceding  Pag^S  Blank 


2 


NACA  TN  I9O3 


Although  these  actions  of  pressure  disturhances  in  aerodynamic 
and  thermodynamic  phenomena  have  become  recognized  and  their  effects 
have  been  partly  evaluated  from  the  general  trend  of  experimental  tests, 
the  amount  of  available  experimental,  data  has  been  insufficient  to 
verify  the  theory  of  unsteady  flow.  This  lack  of  fundamental  experi¬ 
mental  data  has  been  partly  due  to  a  lack  of  instrumentation  capable  of 
responding  with  any  degree  of  accuracy  to  the  rapidly  changing  conditions 
of  this  type  of  flow. 

In  general,  many  of  the  one-dimensional  unsteady— flow  cases  of 
fundamental  interest  can  be  treated  mathematically  with  existing  methods. 
(See  references  3  to  6.)  These  methods,  in  their  present  form,  however, 
are  not  applicable  to  viscous  flows,  flows  involving  heat  transmission, 
or  flows  of  more  than  one  dimension,  which  are  usually  the  flows  of 
engineering  interest. 

The  purpose  of  the  present  investigation  is  to  compare  the  results 
of  an  experimental  study  of  one— dimensional  unsteady  flow  with  results 
based  on  the  theory  of  one-dimensional,  nonviscous,  unsteady  flow.  A 
better  knowledge  of  the  gas  dynamics  of  this  simple  type  of  unsteady 
flow  should  prove  to  be  of  great  assistance  in  the  study  of  the  more 
complicated  flows.  As  a  secondary  objective,  correlation  of  theory  and 
experiment  has  been  used  to  determine  the  value  of  instrumentation  for 
measurements  of  one-dimensional ,  nonviscous,  unsteady  flow. 


SYMBOLS 


p  pressure,  atmospheres 

u  flow  velocity,  feet  per  second 

a  velocity  of  sound,  feet  per  second 

t  time,  seconds 

T  temperature,  degrees  Fahrenheit 

7  ratio  of  specific  heats  at  constant  pressure  to  specific  heats 

at  constant  volume 

x  distance,  inches 

M  Mach  number  in  stationary  shock  coordinates 

V  flow  velocity  in  stationary  shock  coordinates,  feet  per  second 

W  velocity  of  propagation  of  shock,  feet  per  second 
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E  voltage  ,  volts 

l  length.  Inches 

Ap  pressure  change  at  crystal  due  to  given  pressure  front 

Subscripts: 

a  standard  atmospheric  (29.92  in.  Hg) 

0  high-pressure  chamber  before  diaphragm  burst 

1  low-pressure  chamber  before  diaphragm  burst 

b  before  reflection 

c  after  reflection 

22  between  temperature  contact  discontinuity  and  expansion, 

before  reflection 

2  between  temperature  contact  discontinuity  and  shock,  before 

reflection 

3  between  reflected  shock  and  end  of  low-pressure  chamber 

k  between  reflected  expansion  and  end  of  high-pressure  chamber 

L  low-pressure  chamber 

H  high-pressure  chamber 

Superscripts: 

1  high-pressure  side  of  shock 

low-pressure  side  of  shock 

THEORETICAL  METHOD 

The  fundamental  relationships  for  the  propagation  of  plane 
disturbances  of  finite  amplitude  were  first  derived  by  Riemann  in  1859 
for  one-dimensional  isentropic  flow.  (See  references  3  and  7.)  The 
derivation  showed  that  the  propagation  of  a  large  disturbance  can  be 
treated  as  the  propagation  of  an  infinite  number  of  small  disturbances 
where  each  small  disturbance  is  propagated  at  a  velocity  with  respect 
to  the  fluid  equal  to  the  velocity  of  sound  at  that  point.  If  the 
fluid  at  that  point  has  a  flow  velocity  of  u,  then  the  velocity  of 
propagation  of  this  point  on  the  disturbance  is  a  ±  u,  where  a  is 
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the  velocity  of  sound  and  the  plus  sign  or  minus  sign  is  used  if  the 
point  is  "being  propagated  in  the  direction  of  flow  or  against  the 
direction  of  flow,  respectively.  For  a  streamline  of  constant  cross- 
sectional  area  and  for  a  fluid  which  has  a  constant  value  of  7,  a 
disturbance  propagated  in  the  direction  of  flow  was  shown  to  produce 


no  change  in  the  quantity  of  the  flow  u 


through  which  it 


has  been  propagated;  whereas  a  disturbance  propagated  against  the  direc— 

2 

tion  of  flow  was  shown  to  produce  no  change  in  the  quantity  u  +  _  ^  a. 

The  propagation  of  expansion  fronts  on  these  disturbances  was  shown  to 
be  such  as  to  make  the  expansion  fronts  less  steep  with  time  (that  is, 
spreading  expansion  zones);  whereas  compression  fronts  were  shown  to 
become  more  steep  with  time  and  eventually  to  become  discontinuous 
pressure  fronts  (shocks).  From  this  result,  it  can  be  concluded  that 
the  shape  (the  shape  being  the  pressure  distribution  along  the  axis  of 
propagation)  of  any  finite  isentropic  disturbance  continually  changes 
with  time.  In  order  to  calculate  the  shape  of  a  disturbance  at  any 
given  point  in  time,  therefore,  the  shape  must  obviously  be  known  at 
some  earlier  time.  If  an  experimental  measurement  on  the  moving 
disturbance  is  not  trusted,  then  the  time  and  conditions  under  which 
the  disturbance  was  generated  must  be  known  (time  =  0) . 


Classical  methods  for  the  treatment  of  disturbances  from  initial 
conditions  which  are  known  include  the  following: 


1.  Assumption  that  a  piston  rapidly  accelerates  to  a  given  velocity 
which  is  then  held  constant.  This  acceleration  of  the  piston  generates 
a  compression  front  which  moves  ahead  of  the  piston,  while  the  fluid 
behind  the  front  flows  at  the  velocity  of  the  piston.  Deceleration  of 
the  piBton  to  rest  consequently  generates  an  expansion  zone  which  moves 
ahead  of  the  piston  and  reduces  the  fluid  velocity  to  zero. 

2.  Assumption  that  a  diaphragm  which  separates  a  chamber  of  high 
pressure  from  a  chamber  of  low  pressure  is  suddenly  burst.  In  this 
treatment,  a  compression  front  is  propagated  into  the  low-pressure 
chamber;  whereas  an  expansion  zone  is  propagated  into  the  high-pressure 
chamber.  The  movement  of  the  compression  front  is  such  as  to  compress 
the  fluid  in  the  low-pressure  chamber  to  an  intermediate  pressure  and 
to  accelerate  the  compressed  fluid  in  the  direction  of  the  front.  The 
movement  of  the  expansion  zone  is  such  as  to  reduce  the  pressure  in  the 
high-pressure  chamber  to  the  same  intermediate  pressure  and  to  accelerate 
the  expanded  fluid  in  a  direction  opposite  to  the  front.  A  temperature 
contact  discontinuity  appears  behind  the  compression  and  moves  in  the 
direction  of  fluid  flow.  This  surface  of  discontinuity  is  the  boundary 
between  gas  particles  expanded  from  the  high-pressure  chamber  and  particles 
compressed  from  the  low-pressure  chamber. 


The  diaphragm  method  was  chosen  because  it  seemed  somewhat  easier  to 
reproduce  experimentally. 


*- 
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In  1932,  Schardin  (reference  8)  performed  calculations  from  which 
the  pressure  distribution  along  tubes  of  constant  cross  section  can  be 
found  at  certain  times  after  the  shattering  of  the  diaphragm.  These 
calculations  are  based  on  the  assumptions  that  the  flow  is  nonviscous 
and  that  a  plane  shock  wave  is  immediately  formed  which  moves  into  the 
tube  of  low  pressure  while  an  expansion  wave  moves  into  the  tube  of  high 
pressure.  Schardin  used  the  Eankine-Eugoniot  relations  for  the  shock 
fronts  and  Riemannfs  relations  for  the  expansion  fronts. 

By  an  extension  of  Schardin*  s  calculations  to  include  reflections 
and  interactions  of  the  waves,  as  well  as  the  effect  of  the  temperature 
discontinuity  which  exists  in  the  flow,  a  complete  time  history  of  the 
phenomena  can  be  calculated.  In  order  to  find  the  conditions  which  exist 
immediately  after  the  shattering  of  the  diaphragm,  the  method  used  by 
Schardin  should  be  briefly  reviewed. 

In  figure  1  is  shown  the  pressure  distribution  at  typical  times 
along  constant— area  passages  which  were  separated  by  "the  diaphragm  at 
t  =  0.  Figure  l(a)  shows  the  pressure  distribution  along  a  tube  in 
which  a  diaphragm  separates  a  pressure  Pq  from  a  lower  pressure  p^« 

Figure  1(b)  shows  the  pressure  distribution  along  the  tube  at  a  time  t^ 

after  the  diaphragm  has  been  shattered.  A  shock  front  is  assumed  to  be 
propagated  towards  the  region  of  low  pressure  and  to  be  compressing  the 
gas  at  p^  to  an  intermediate  pressure  pg.  At  the  same  time,  a  zone 

of  expansion  is  assumed  to  be  propagated  towards  the  direction  of  high 
pressure  and  to  be  reducing  the  pressure  p q  to  the  same  intermediate 

pressure  p^.  In  order  to  find  all  the  flow  quantities  existing  in  the 
tube  at  the  time  t-^,  the  following  assumptions  are  made: 

(1)  The  value  of  7  is  constant. 

(2)  Velocities  are  positive  in  the  positive  x-direetion. 

(3)  The  tube  is  of  constant  cross-sectional  area. 

(4)  Before  the  diaphragm  is  shattered,  the  temperature  in  the 

tube  is  everywhere  the  same  and  the  velocity  is  everywhere 
equal  to  0. 

(5)  The  flow  is  nonviscous  and  involves  no  transmission  of  heat. 

Immediately  after  the  rupturing  of  the  diaphragm,  a  temperature 
contact  discontinuity  exists  in  the  flow  somewhere  between  the  compression 
front  and  the  expansion  front  as  mentioned  previously.  This  discontinuity 
also  represents  an  entropy  discontinuity,  because  this  discontinuity 
represents  no  pressure  increment  and,  hence,  no  velocity  increment.  In 
order  to  satisfy  the  assumption  that  an  expansion  zone  is  moving  in  the 


r 
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negative  x-direction,  it  now  follows  from  Riemann’s  work  that  the 
velocity  U2  must  he  in  the  positive  x— direction.  Therefore, 

r 


U  0  + 


Ug  + 


a22 


which  may  be  expressed  as 


u2  -  u0 
a0 


(1) 


because  the  process  is  isentropic  and 

7-1 

*22  _  /PA27 
a0  W 

Equation  (l)  expresses  the  velocity  increment  across  the  expansion  zone 
as  a  function  of  the  pressure  ratio  across  the  zone  and  aQ,  which  is  a 
function  of  the  initial  temperature. 

The  velocity  increment  across  the  shock  front  can  be  expressed  in 
terms  of  the  pressure  ratio  across  the  front  and  a^  by  using  the 

Rankine-Hugoniot  relations  for  a  stationary  normal  shock  as  follows. 

The  coordinate  system  in  figure  l(b)  can  be  changed  to  reduce  the 
shock  velocity  =  0  by  expressing  the  gas  velocities  on  either  side 

of  the  shock  as 


V" 


Un  -  W, 


* 


QJll 


f 


V  =  u2  -  Wb 
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where  V"  is  the  gas  velocity  on  the  low-pressure  side  of  a  stationary 
normal  shock  and  V*  is  that  on  the  high— pressure  side.  (See  fig.  1(c).) 
The  velocity  increment  across  the  shock  is  then 


V*  -  Y"  =  -M2a2  -  (-%&!) 


which  may  he  written  as 


Y»  -V" 
al 


(2) 


where  and  M2  are  M8*0*1  numbers  on  the  low-pressure  and  high- 
pressure  side,  respectively,  of  a  stationary  normal  shock  and  from  the 
Rankine-Bugoniot  relations. 


Also,  from  the  stationary— shock  relations. 


%  = 


(7  + 1)—  +  (7-1) 

pi 


27 


and 


'(7  +  l)^  +  (7-1) 

p2 


M g  = 


27 
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* 


so  that  equation  (2)  may  be  -written  as 


Equation  (3)  is  the  expression  for  the  velocity  increment  across  the 
shock  front,  for  either  coordinate  system,  as  a  function  of  the  pressure 
ratio  across  the  front  and  ali  It  might  be  pointed  out  that  the  selection 
of  coordinate  systems  other  than  the  one  used  herein  yields  the  same 
relation.  (For  example,  see  reference  5.)  Equation  (3)  may  be  equated 
to  equation  (l)  because  aq  =  aQ  ^d  U1  =  u0  =  0  as  follows: 


Pi  =  Pi 

p0  P2 


-|2  7 


7-1 


ft 


w 


Equation  (4)  can  easily  be  solved  graphically  for  P2/Pq.  (See 
fig.  2.)  By  substituting  p2/p^  in  equation  (3),  Ug  may  be  found 
because  aq  is  known  frcrn  the  initial  temperature  in  the  tube.  Also, 
Wq  may  be  found  from 


-  (“Mlal) 


or  the  pressure  ratio  across  the  shock  may  be  written  as  a  function  of 
the  shock  speed 


=—  I  —  1 


7  +  1 


(5) 


« 

/ 


when  uq  =  0. 
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All  the  conditions  in  the  tube  are  now  known  for  any  time  t-^  after 
the  shattering  of  the  diaphragm  -until  reflections  occur  at  the  ends  of 
the  tube.  These  conditions  satisfy  the  original  assumptions  of  an 
expansion  zone  moving  in  the  negative  x— direction  and  a  compression  shock 
moving  in  the  positive  x-direction,  aB  well  as  a  temperature  contact 
discontinuity  (a22  i  moving  a*k  the  velocity  u2  across  which  there 
is  no  change  in  pressure  or  flow  veloqity.  The  pressure  distribution 
along  the  tube  is  easily  found  because  the  shock  has  moved  a  distance.  W^t^ 
and  the  expansion  zone  AB  in  figure  l(b)  has  moved  distances  of  -e-o^b 
and  —  (a2p  —  u2^b  ^or  points  A  and  B,  respectively.  The  expansion 

zone  becomes  less  steep  with  time  because  point  A  moves  faster  than 
point  B;  whereas  a  compression  shock  exists  immediately  and  moves  in  the 
positive  x-direction  because  the  propagation  velocity  behind  the 
shock  (u£  +  a^)  is  greater  than  that  ahead  (ap). 

Figure  l(d)  shows  the  pressure  distribution  along  the  tube  at  a 
time  tc  after  the  shock  has  been  reflected  from  the  low-pressure  end  of 
the  tube  and  the  expansion  has  been  reflected  from  the  other  end.  The 
conditions  existing  after  reflection  of  the  shock  can  be  found  in  a 
manner  similar  to  that  used  to  obtain  equation  (3);  that  is,  the  reflected 
shock  velocity  — Wc  =  0  can  be  reduced  so  that  the  flow  velocities  on 
either  side  of  the  shock  become  (fig.  1(e)) 


v”  =  u2  -  (_wc) 

and 

V»  =  u3  -  (-tfc) 


The  velocity  increment  across  the  shock  can  therefore  be  expressed  in 
terms  of  its  pressure  ratio  and  ^  a8  follows: 


yt  _  y,f  u^  —  U2 

a2  a2 
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Because  U3  =  0  at  the  end  of  the  tube,  P3/P2  can  expressed 
as  a  function  of  Pg/Pq  by  combining  equation  (6)  with  equation  (3)  where 


The  resulting  expression  is,  for  7  =  l-.U, 


2 


2 


h 

p2 


1&-1 

?! 

^  +  6 

pi 


(7) 


The  shock  Telocity  Wc  can  be  found  from 


-¥c  =  u3  -  V* 


or 

-*c  =  *2  -  V"  =  V2 


where  and  are  the  Mach  numbers  on  the  high-pressure  and  low- 

pressure  side,  respectively,  of  the  reflected  shock  in  stationary 
coordinates;  that  is. 


Mn  = 


'7  +  1  £2 

=Ah 

27... 


+ 1 


7-1 


= 


[7  + 1 P3 


7  -  1  P2 
27 


+  1 
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and 


f3 

a2 


All  the  flov  conditions  after  the  reflection  of  the  shock  are  now 
known.  Schardin’s  calculations  (reference  8)  are  extended  to  the  reflec¬ 
tion  of  the  shock.  (See  table  I.)  , 

In  order  to  extend  the  calculations  to  a  complete  history  of  the 
flow  inside'  the  tube,  the  next  step  is  to  find  the  reflection  of  the 
expansion  zone  from  the  end  of  the  tube.  The  conditions  after  total 
reflection  of  the  expansion  are  easily  found  from  Eiemann’s  equation 
where 

u4  -  Y&J  a4  =  u2  -  yzrj  a22 

which  may  be  solved  for  a^  because  u^  =  0  at  the  end  of  the  tube. 

The  pressure  p^  can  then  be  found  from  the  isentropic  relationship 

27 

P4  =  /v\7-1 

P2  ~  \a22/ 

Because  of  its  thickness,  however,  the  total1  reflection  of  the 
expansion  requires  a  finite  time  so  that  the  pressure  history  at  the 
end  of  the  tube  during  reflection,  or  the  pressure  distribution  along 
the  axis  of  propagation  between  points  A  and  B  must  be  found  somewhat 
differently.  This  pressure  history  is  obtained  by  considering  the 
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expansion  zone  AB  to  consist  of  an  arbitrary  number  of  small  discontinuous 
isentropic  expansion  fronts  (the  usual  characteristics  method)  rather 
than  to  be  a  smooth  pressure  variation  between  zone  AB,  For  an  illustra¬ 
tion  of  this  method,  see  figure  3  in  which  the  expansion  zone  AB  is 
assumed  to  consist  of  two  expansion  zones  Y  and  Z  which  are  replaced  with 
isentropic  fronts  Y  and  Z.  This  method  is,  of  course,  not  rigorous  because 
Eiemann  has  shown  that  an  expansion  front  cannot  remain  discontinuous; 
however,  by  the  choice  of  small  Bteps  the  method  provides  a  good  approxima¬ 
tion,  Each  front  can  now  be  considered  individually,  and  the  conditions 
of  u  and  a  on  either  side  of  each  front  can  be  found  from  the  initial 
pressure  increment  assigned  the  front  by 


u 


+ 


a  =  Constant 


In  figure  3,  step  2,  p  =  1.70  atmospheres  was  arbitrarily  assigned 
the  region  between  fronts  Y  and  Z.  From  the  relation  u  +  5a  =  Constant 
across  either  front  Y  or  Z  (because  the  fronts  are  being  propagated 
against  the  fluid  flow)  and  the  adiabatic  relations,  the  value  of  u 
and  a  in  this  region  are  found.  The  speed  of  each  front  is  assumed  to  # 

be  an  average  of  the  values  of  u  ±.  a  on  either  side  of  the  front.  For 
this  reason,  the  characteristics  method  represents  an  approximation. 

The  total  reflection  of  the  expansion  from  the  end  of  the  tube  is 
obtained  by  reflection  of  the  small  discontinuous  expansion  fronts 
individually  and  by  consideration  of  the  subsequent  interactions  of  the 
reflected  with  unreflected  fronts.  In  figure  3,  step  3,  the  conditions 
at  the  tube  end  after  reflection  of  front  Y  are  first  found  by  using  the 
relation  u  —  5 a  across  reflected  front  Y  and  by  knowing  that  u  =  0  at 
the  tube  end.  The  conditions  between  fronts  Y  and  Z  after  interaction 
are  next  found  by  using  the  relationships  u  +  5&  =  Constant  across 
front  Z  and  u  —  5a  =  Constant  across  front  Y  and  by  obtaining  a 
simultaneous  solution  for  u  and  a  between  fronts  after  the  interaction. 

(See  fig.  3,  step  4.)  The  final  step  is  the  reflection  of  front  Z 
from  the  tube  end.  (See  fig.  3j  step  5-) 

Because  the  speeds  of  all  the  fronts  are  known  and  because  all  the 
fronts  originate  at  the  diaphragm  at  t  =  0,  a  more  convenient  method  of 
illustrating  the  positions  of  the  various  fronts  in  the  tube  at  any  time 
is  to  plot  the  time  histories  as  a  diagram  of  x  against  t  as  in 
figure  4  where  the  expansion  is  divided  into  seven  fronts.  Pressure 
variations,  as  a  function  of  time  or  distance  along  the  tube,  are  easily 
taken  from  such  a  plot  because  conditions  in  each  enclosed  area  are 

known  and  are  constant  throughout  the  area.  In  figure  5,  the  preceding  * 

method  was  used  to  show  the  pressure  distribution  along  the  tube  between 

points  A  and  B  at  a  time  t^.  In  figure  6,  this  method  was  used  to  show 

the  variation  of  pressure  at  the  end  of  the  tube  with  respect  to  time  at  / 

reflection  of  expansion  zone.  It  can  be  seen  from  figures  5  ancl  6  that 

the  pressure  variation  inside  the  expansion  front  is  nonlinear. 
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Figure  4  shows  that,  depending  on  the  ratio  of  the  lengths 

%  the  temperature  contact  discontinuity  is  overtaken  hy  either  the  reflected 

expansion  zone  or  the  reflected  shock  front.  Whichever  is  the  case,  it 
can  he  shown  that  the  pressure  ratio  of  the  front  is  altered  after  passing 
through  the  temperature  contact  discontinuity  which,  in  turn,  generates 
a  second  front  that  propagates  itself  in  the  opposite  direction.  The 
method  used  to  obtain  the  conditions  which  exist  after  the  front  has 
undergone  this  interaction  with  the  temperature  contact  discontinuity  is 
similar  to  the  method  already  used.  First,  the  assumption  is  made  that 
a  second  front  is  generated  and  travels  in  the  opposite  direction  from 
the  original  front.  Then,  knowing  that  there  is  no  velocity  increment 
or  pressure  increment  across  the  new  discontinuity  necessitates  the 
finding  of  four  unknown  conditions:  namely,  the  pressure  and  velocity 
existing  in  the  region  between  the  two  fronts,  and  the  two  temperatures. 

These  conditions  can  be  found  in  the  same  manner  in  which  the  condi¬ 
tions  after  the  shattering  of  the  diaphragm  were  found,  except  that  the 
generated  front  must  be  assumed  to  be  either  a  shock  or  an  expansion. 

After  solutions  are  obtained  and  the  assumption  is  indicated  to  be 
^  incorrect,  a  second  solution  must  be  obtained.  If  the  interaction  under 

consideration  is  that  of  an  expansion  zone  with  a  discontinuity,  the 
expansion  zone  must  be  considered  as  a  number  of  small  fronts  (as  was 
the  treatment  for  the  reflection  of  an  expansion  zone  from  the  end  of 
the  tube)  in  order  to  obtain  a  pressure  variation  with  time  or  distance 
inside  the  expansion  zone  during  the  interaction.  This  method  is  rather 
laborious  and  is  not  required  when  solving  for  conditions  after  the 
interaction  but  should  be  used  for  an  accurate  history  of  events  in  a 
diagram  of  x  against  t.  Interactions  of  expansion  zones  and  shock 
fronts  with  discontinuities  have  been  analyzed  by  this  method,  and,  in 
general,  the  generated  front  is  very  weak  in  proportion  to  the  original 
front.  An  interaction  of  a  shock  front  with  a  discontinuity  is  shown  in 
figure  4  in  which  the  generated  front  is  a  compression  traveling  in  the 
opposite  direction  from  the  original  shock. 

In  order  to  be  able  to  find  conditions  in  the  tube  at  any  time  or 
position,  the  one  remaining  phenomenon  which  must  be  considered  is  an 
interaction  of  a  shock  front  with  an  expansion  zone.  For  this  case, 
the  conditions  which  exist  after  the  interaction  is  over  are  again 
readily  found  by  the  methods  used  for  finding  the  conditions  after  the 
shattering  of  the  diaphragm.  The  only  assumption  that  satisfies  a 
solution  of  this  interaction  is  one  in  which  a  shock  and  an  expansion 
are  propagated  in  the  same  direction  after  the  interaction  as  before 
the  interaction  and  in  which  the  flow  between  the  two  fronts  is  reversed 
in  direction.  The  values  of  u,  a,  and  p,  which  exist  between  the 
fronts  after  interaction,  can  be  obtained  by  making  this  assumption  and 
by  obtaining  simultaneous  solutions  of  the  shock  equations  and  the 
isentropic— front  equations.  The  expansion,  however,  must  be  considered 
'  as  a  number  of  small  fronts  and  the  method  must  be  applied  the  same 

number  of  times  if  a  pressure  variation  of  time  or  distance  is  to  be 
obtained  during  the  interaction.  An  interaction  between  a  shock  front 
and  an  expansion  zone  is  shown  in -figure  4. 
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By  use  of  the  methods  which  have  heen  discussed  it  is  possible  to 
determine  all  the  subsequent  phenomena  and  to  obtain  the  complete  pressure 
history  in  the  tube  at  any  point  in  the  tube  or  a  complete  pressure 
distribution  along  the  tube  at  any  time  after  the  diaphragm  is  shattered. 
The  assumptions  of  nonviscous  flow  of  a  gas  which  has  a  constant  value 
of  7  and  which  involves  no  heat  transmission  should  be  valid  in  many 
cases.  However,  for  passages  of  very  small  cross  section  or  long  length 
and  for  pressure  variations  involving  large  temperature  variations, 
these  assumptions  would  undoubtedly  lead  to  incorrect  results.  The 
method  which  has  been  presented  can,  in  many  cases,  be  greatly  simplified 
for  small  pressure  variations.  In  these  cases  the  shock  can  be  treated 
as  an  isentropic  pressure  front,  and  Riemann* s  relations,  rather  than  the 
shock  relations,  can  be  used  with  very  little  error.  Also,  the  expansion 
zone  need  not  be  broken  into  a  large  number  of  weaker  fronts;  usually  two 
or  three  will  suffice. 


EXPERIMENTAL  METHOD 


A  diagrammatic  sketch  and  photographs  of  the  apparatus  are  shown  in 
figures  7  and  8,  respectively.  A  description  of  the  design  and  function 
of  the  instrumentation,  much  of  which  was  developed  specifically  for 
these  tests  by  the  Langley  Instrument  Research  Division,  was  prepared 
by  Mr.  E.  Delpino  and  is  given  in  the  appendix. 

The  apparatus  in  which  the  waves  were  produced  consisted  of  four 
passages  of  constant  Internal  cross-sectional  area  which  could  be 
connected  together  in  any  sequence  by  means  of  clamps  and  by  mating 
flanges  to  form  one  long  passage  with  smooth  concentric  inner  surfaces. 
The  internal  dimensions  were  everywhere  3  inches  wide  by  3  inches  high. 

The  walls  were  fabricated  from  ^--inch  brass  stock  by  brazing  the  joints 

along  the  outside  surface  of  the  passage.  The  four  sections  included  a 
section  8  inches  long  which  contained  a  solenoid— actuated  trigger  used 
to  puncture  the  diaphragm.  (See  fig.  7(c).)  Also  included  was  a  2^ — inch 

section,  a  48-inch  section,  and  a  24“  inch  section,  which  was  equipped 

with  two  pairs  of  glass  windows  of  such  dimensions  that  the  full  3-inch 
height  of  the  passage  for  a  distance  of  6  inches  along  the  passage 
could  be  photographed  from  either  pair  of  windows.  The  glass  windows 
were  mounted  in  such  a  manner  that  the  inner  surfaces  of  the  glass  were 
flush  with  the  inner  surfaces  of  the  passage  with  gaskets  between  the 
joints  to  prevent  leakage  of  air.  Each  section  of  passage  was  equipped 
with  two  small  static— pressure  orifices  (0.04— inch  diameter),  one  which 
could  be  connected  to  a  pressure  gage  and  the  other,  to  a  source  of 
pressure  or  vacuum.  At  the  positions  shown  in  figure  7(n),  each  section 
was  equipped  with  tapped  holes  along  the  top  of  the  passage  which  could 
accomodate  either  a  blank  plug  screwed  in  flush  with  the  inner  surface 
or  a  plug  in  which  the  pressure-measuring  instrument  was  mounted  flush 
with  the  surface  as  shown  in  figure  7(b).  Three  metal  end  plates  were 
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provided  of  thicknesses  of  3/16  inch,  3/8  inch,  and  11 /l6  inch,  each  of 
which  could  he  clamped  to  the  end  of  any  section  of  passage. 

In  all  the  tests,  the  8— inch  trigger  section,  the  2*1— inch  section, 

and  the  5-  — inch  end  plate  were  clamped  together  to  form  what  is  called 
o 

the  high-pressure  chamber.  The  low-pressure  chamber  then  consisted  of 

either  or  both  of  the  remaining  two  sections  in  addition  to  the  11.—  inch 

l6 

end  plate.  In  this  manner  three  different  lengths  of  low-pressure  chamber 
were  available  in  four  possible  combinations  of  section  —  that  is,  one 

2^1  inch,  one  k8  inch,  and  two  72^-  inch  —  in  which  the  glass  section 

could  be  at  either  end.  All  the  clamped  Joints  between  sections  and 
between  a  section  and  an  end  plate  were  equipped  with  rubber  gaskets  of 
such  dimensions  that  the  inner  edges  were  flush  with  the  inner  surfaces 
of  the  passages.  The  trigger  was  equipped  with  a  sliding  pressure  seal 
and  a  soft  rubber  cushion  to  damp  out  mechanical  vibration.  (See  fig.  7«) 
The  blank  plugs  and  instrument  plugs  were  equipped  with  rubber  gaskets. 

The  diaphragm  used  to  separate  the  high-pressure  chamber  from  the  low- 
pressure  chamber  consisted  of  from  one  to  four  thicknesses  of  a  special 
paper,  depending  on  the  pressure  difference  between  the  two  chambers, 
and  was  inserted  in  this  Joint  along  with  a  rubber  gasket  before  clamping. 
The  particular  paper  used,  American  Blueprint  Co.  Inc.,  No.  *J-5D,  Type  2, 
tracing  paper,  was  selected  after  many  preliminary  trials  because  of  the 
following  favorable  characteristics:  very  little  stretching  under 
tension;  minimum  porosity  and  hence  minimum  leakage  of  air  through  the 
paper  before  puncturing;  and,  finally,  disintegration  into  smal  1  pieces 
when  punctured  rather  than  simply  tearing  into  four  pie-shaped  pieces. 

The  pressure-measuring  instrument  was  a  commercial  unit  consisting  of 
an  ammonium  dihydrogen  phosphate  piezoelectric  crystal  mounted  in  a  metal 

case  of  S.  — inch  diameter  and  1— inch  length.  The  electrical  characteristics 
o 

of  this  instrument  are  described  in  the  appendix.  This  instrument  has 
the  following  desirable  characteristics  for  this  type  of  test:  rugged 
construction,  high  acoustic  impedance,  high  sensitivity,  and  no  cavity 
resonance. 

Before  each  series  of  tests,  the  high-pressure  and  low-pressure 
chambers  were  clamped  together  without  a  diaphragm  and  the  entire  unit 
was  checked  for  air  leakage.  Before  each  test,  the  small  bits  of  paper 
left  in  the  chambers  from  the  preceding  test  were  blown  out  with  an  air 
hose.  After  a  new  diaphragm  was  clamped  between  the  high— pressure  and 
low-pressure  chambers,  the  pressure  in  each  chamber  was  adjusted  to  the 
desired  values  and  read  from  calibrated  pressure  gages  to  an  accuracy 
of  0.1  pound  per  square  inch.  The  pressures  wete  held  Bteady  for  a  few 
minutes  before  each  test  so  that  the  temperatures  could  come  into 
equilibrium  with  the  room  temperature.  Because  the  gage  leads  and 
pressure— supply  leads  were  connected  to  the  apparatus  through  small 
orifices  and  because  the  duration  of  the  test  was  only  a  few  milliseconds. 
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the  small  amount  of  flow  Induced  through  these  orifices  during  the  test 
were  "believed  to  "be  of  no  consequence. 

Photographs  of  the  flow  were  taken  with  a  conventional  schlieren 
system  which  included  two  parabolic  front— surfaced  mirrors  of  6-inch 
diameter  and  a  high— voltage— spark  light  source  which  was  synchronized 
with  the  diaphragm  hurst  as  described  in  the  appendix#  Measurements  of 
the  absolute  magnitude  of  the  shock  speed  before  its  reflection  from 
the  end  of  the  low-pressure  chamber  were  made  as  described  in  the  appendix# 

Preliminary  tests  were  made  by  using  a  setup  of  smaller  dimensions 
than  the  one  previously  described.  These  tests  were  made  principally  in 
order  to  develop  the  instrumentation.  Additional  preliminary  tests  were 
made  by  using  the  apparatus  described  in  order  to  determine  the  effects 
of  the  following  objectionable  conditions  on  the  pressure  records: 

(1)  Air  leakage  or  air  cavity  around  the  crystal' unit 

(2)  Mounting  of  the  crystal  unit  in  rubber  instead  of  the  Micarta 
mount  shown  in  figure  7(b) 

(3)  Mechanical  vibration  caused  by  the  solenoid-actuated  trigger 

(4)  Failure  of  the  crystal  unit  to  be  perfectly  flush  with  the 
inner  surface  of  the  passage  or  end  plate 

(5)  Failure  of  the  diaphragm  to  disintegrate  after  bursting 

(6)  Mechanical  vibration  of  the  end  plates  set  up  by  the  reflection 
of  the  shock  waves 

Before  the  effect  of  these  conditions  on  the  test  results  are 
discussed,  the  method  of  interpreting  the  test  records  should  be  discussed# 
As  pointed  out  in  the  appendix,  the  test  records  are  essentially  a  plot 
of  the  voltage  developed  by  the  crystal  unit  as  a  function  of  time 
provided  that  proper  electrical  instrumentation  is  used.  If,  therefore, 
the  pressure  sensitivity  of  the  crystal  unit  is  known  (that  is.  the 
voltage  output  as  a  function  of  pressure  change  at  the  crystal),  the 
test  results  can  then  be  compared  with  theory  by  comparing  the  theoreti¬ 
cal  and  experimental  pressure— time  plots.  The  theoretical  pressure— time 
plot  for  any  position  of  the  crystal  unit  can  be  taken  from  plots  such 
as  those  shown  in  figures  9(a)  to  9(c),  which  were  obtained  for  a  pressure 
ratio  Pq/Pi  2.0  for  the  three  lengths  of  low-pressure  chamber,  by 

using  the  methods  outlined  previously  and  by  assuming  the  fluid  to  be 
air  (7  =  1.4).  In  figure  10,  the  theoretical  pressure— time  plots  are 
shown  as  dash  lines  for  a  number  of  different  locations  and  lengths  of 
low-pressure  chamber,  all  for  Pq/Pi  °f  2.0. 

With  the  theoretical  pressure— ti me  plot  known,  the  effects  of  the 
objectionable  conditions  previously  listed  were  experimentally  studied 
and  the  results  of  these  preliminary  tests  are  as  follows: 
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(1)  Records  were  taken  with  and  without  air  leakage  and  an  air 
cavity  around  the  crystal  unit.  Comparison  of  the  records  indicated  that 
a  sharp  fall  of  the  trace  occurred  Immediately  following  a  pressure  rise 
at  the  crystal  when  leakage  and  air  cavity  were  present. 

(2)  Records  were  taken  with  the  crystal  unit  mounted  in  rubber  and 
were  compared  with  records  taken  with  the  unit  mounted  in  the  Micarta 
mount  shown  previously.  Oscillations  of  the  trace  on  the  records  taken 
with  the  rubber  mount  occurred  following  a  sharp  compression.  These 
oscillations  indicated  that  the  crystal  was  oscillating  mechanically  at 
the  frequency  of  the  rubber  mount. 

(3)  Records  were  taken  with  no  pressure  difference  across  the  diaphragm, 
and  a  comparison  was  made  of  the  trace  obtained  with  and  without  a  soft 
rubber  cushion  on  the  trigger  mechanism  shown  previously.  The  mechanical 
vibration  caused  by  the  trigger  hitting  the  stop  on  the  side  of  the 
passage  was  found  to  be  substantially  eliminated  by  the  rubber  cushion. 

(k)  Records  taken  with  the  surface  of  the  crystal  unit  not  flush  with 
the  inner  surface  of  the  passage  or  end  plate  indicated  that  the  trace 
was  distinctly  altered  because  of  the  waves  reflected  from  the  crystal 
edges  in  this  condition. 

(5)  Records  were  taken  of  bursts  in  which  more  thicknesses  of  paper 
diaphragm  were  used  than  necessary,  in  which  case  the  burst  was  more  a 
tearing  of  the  paper  than  a  disintegration  of  the  paper  into  tiny  pieces. 
These  results  indicated  that  in  cases  where  the  paper  did  not  disintegrate 
the  rise  of  the  trace  as  a  result  of  a  compression  was  more  gradual  and 
was  followed  by  oscillations  of  the  trace,  probably  a  result  of  small 
waves  behind  the  compression  front  caused  by  poor  air  flow  past  the 
remaining  part  of  the  diaphragm. 

(6)  On  all  records  taken  with  the  crystal  unit  in  the  end  plates,  a 
rather  violent  oscillation  of  the  trace  would  always  occur  following  a 
compression.  This  oscillation  would  not  occur  follbwing  an  expansion, 
nor  would  it  occur  for  either  case  if  the  crystal  unit  was  mounted  in  the 
side  of  the  passage.  Also,  the  use  of  thicker  end  plates  increased  the 
frequency  of  this  oscillation  but  had  little  effect  on  its  amplitude. 

The  addition  of  a  ^-pound  weight  onto  the  plug  in  which  the  crystal  was 
mounted  (fig.  7(h))  would  decrease  the  frequency  of  oscillation  for.  a 
given  end— plate  thickness.  In  all  cases  the  rise  of  the  trace  due  to 
the  compression  was  reproduced  (for  a  given  pressure)  if  measured  to  the 
mean  of  this  oscillation.  The  oscillation  was  a  result  of  mechanical 
vibration  of  the  end  plate  excited  by  the  sharp  pressure  change  in  the 
compression  and  was  extremely  difficult  to  eliminate. 

The  qualitative  nature  of  all  these  effects  can  be  predicted  from 
general  unsteady— flow  theory;  however,  it  was  desired  to  determine  the 
order  of  magnitude  of  the  effects  before  analyzing'the  records.  As  a 
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result  of  these  prellmliiary  tests,  the  standard  conf  iguration  for  all 
subsequent  tests,  unless  otherwise  noted,  incorporated  a  Micarta  crystal 
unit  mount  as  shown  in  figure  7(h),  a  soft  rubber  cushion  on  the  trigger 

as  shown  in  figure  7(c),  an  ^  — inch— thick  end  plate  on  the  low-pressure 

chamber,  and  —  as  pointed  out  in  the  appendix  —  a  100-micromicrofarad 
condenser  across  the  crystal  as  well  as  a  ^7-kilocycle  filter  in  the 
preamplifier  circuit. 

RESULTS  AND  DISCUSSION 


A  large  number  of  test  records  were  taken  at  various  crystal  positions 
and  with  the  three  different  low— pressure— chamber  lengths.  Records  were 
taken  over  a  range  of  pressure  ratio  Pq/pi  across  the  diaphragm  from  1.3 

to  2.5  where  p-^,  the  pressure  in  the  low-pressure  chamber,  was  equal 

to  1  atmosphere  or  l/2  atmosphere.  Topical  pressure  records  which  were 
obtained  are  shown  in  figure  11. 

Before  the  experimental  records  can  be  superimposed  on  a  given 
pressure— time  scale  for  comparison  with  theoiy,  it  is  necessary  to  know 
the  voltage  output  of  the  crystal  unit  as  a  function  of  the  pressure 
change  experienced  by  the  crystal  unit.  This  pressure  sensitivity  of 
the  crystal  unit  was  found  experimentally  by  plotting  all  the  voltage 
changes  as  measured  from  the  test  records  as  a  function  of  the  pressure 
changes  experienced  by  the  crystal  as  predicted  by  the  theories  previously 
outlined.  In  figure  12  and  table  II  are  shown  the  results  of  these 
voltage  measurements,  along  with  the  positions  on  the  test  records  at 
which  voltage  measurements  were  taken.  These  data  indicate  that,  with 
only  a  very  small  scatter  of  points,  a  straight  line  can  be  drawn  through 
the  points  starting  from  zero;  that  is,  the  output  of  the  crystal  unit 
is  a  linear  function  of  the  pressure  change  at  the  crystal,  within  the 
range  of  pressures  shown.  As  shown  in  figure  13*  the  voltage  measure¬ 
ments  E  were  taken  from  the  test  records  by  first  drawing  a  straight 
line  (shown  as  a  dash  line)  through  the  part  of  the  trace  immediately 
following  either  an  expansion  or  compression  and  extending  it  to  the  point 
in  time  at  which  the  voltage  change  began.  This  line  corrects  the 
measurement  for  the  slight  fall  of  the  trace  towards  zero  voltage  inherent 
in  the  crystal  and  electrical  circuit  as  pointed  out  in  the  appendix. 

No  significant  difference  was  noted  among  records  taken  with  a  number  of 
crystals  all  of  the  same  model  type. 

Because  the  abcissa  in  figure  12  is  the  theoretical  pressure  * 

difference  experienced  by  the  crystal  unit  for  a  given  P0/Pi  across 

the  diaphragm,  it  was  necessary  to  determine  experimentally  if  the 

absolute  magnitude  of  this  pressure  difference  was  that  predicted  by  ' 

theory.  In  order  to  do  this  without  making  a  pressure  measurement,  the 
value  of  the  shock  speed  was  experimentally  measured  (see  appendix) 

and  the  pressure  pg  was  calculated  from  equation  (5)*  Pq  al 
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"being  known.  The  results  of  these  measurements  are  shown  in  figure  2, 
where  it  can  be  seen  from  the  measured  values  of  Pg/Pi  that  the  actual 

values  of  Ap  are  only  a  few  percent  lower  than  the  theoretical  values. 
The  use  of  Ap  to  calculate  the  shock  speeds  gave  values  only  a  fraction 
of  a  percent  lower  than  theoretical  shock  speeds. 


p0 

The  experimental  pressure— time  records  for  ^=2.0  were  supers 

imposed  on  the  theoretical  pressure— time  plots  and  are  shown  as  solid  lines 
in  figure  10.  In  figure  10  the  pressure  scales  (that  is,  voltage  scales) 
for  the  theoretical  trace  were  taken  from  figure  12,  and  the  test  records 
were  superimposed  by  tracing  directly  from  the  records  on  a  light  table. 

The  time  scales  and  voltage  against  trace  deflection  scales  were  taken 
from  the  timing  frequency  traces  and  the  alternating— current  voltage  traces 
on  the  test  records. 


Figure  10  shows  that  the  agreement  between  the  theoretical  and 
experimental  pressure— time  plots  is  very  close.  The  slight  departure 
which  is  shown  is  seen  to  be  mainly  due  to  the  fall  of  the  trace  towards 
zero  voltage,  previously  mentioned  as  a  condition  inherent  in  the  electrical 

circuit.  This  fall  was  measured  (see  appendix)  and  was1  about  l|-  percent 
of  the  voltage  per  millisecond.  Although  this  correction  was  not  made  in 
figure  10,  the  agreement  would  be  still  closer  if  the  correction  were  made. 

Schlieren  photographs  taken  by  using  the  2^  — inch  low-pressure 
chamber  and  by  using  the  72g-  —inch  low-pressure  chamber  with  the  glass 

section  at  the  end  plate  are  shown  in  figure  14.  As  a  result  of  the 
schlieren  photographs,  shock— speed  measurements,  and  pressure  records,  a 
number  of  interesting  comparisons  can  be  made  between  the  experimental 
results  and  the  theoretical  predictions  which  were  based  on  one-dimensional 
unsteady  flow. 

The  first  assumptions  made  were  that  the  flow  was  frictionless  and 
that  a  plane  compression  shock  is  immediately  formed  after  the  shattering 
of  the  diaphragm.  In  figure  12,  points  taken  at  crystal  locations 

ranging  from  7?y  inches  from  the  diaphragm  up  to  72“  inches  from  the 

diaphragm  indicate  that  the  shock  does  not  exhibit  any  measurable 
attenuation  within  the  scatter  of  the  points  (less  than  1  percent  of  Ap). 
This  assumption  is  further  confirmed  by  the  pressure— time  records  in 
figure  10  and  the  shock— speed  measurements  in  figure  2  where  attenuation 
would  show  up  as  a  change  in  shock  speed  and  thereby  cause  disagreement 
with  the  theoretical  pressure-time  plot.  The  schlieren  photographs  and 
pressure— time  records  show  that  the  compression  is  not  immediately  one 
dimensional  and  discontinuous  but  that  it  becomes  so  very  rapidly. 

At  inches  from  the  diaphragm  the  schlieren  photographs  and  pressure- 

time  records  (figs.  l4(a),  l4(b),  and  10(d))  indicate  that  the  compres¬ 
sion  has  a  definite  thickness;  however,  the  schlieren  photographs  and 
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pressure— time  records  taken  at  17  inches  and  24^-  inches,  respectively, 

(figs.  l4(d),  10(a),  10(b),  and  10(c))  indicate  that  the  shock  is  very 
sharp  after  this  length  of  travel.  This  rapid  steepening  is  to  be 
expected  from  theoretical  considerations  because  for  this  Pq/p^  the 

value  of  u  +  a  behind  the  shock  is  appreciably  larger  than  that  ahead 
of  the  shock  (l462  fps  as  compared  with  1126  fps).  It  cannot  be  definitely 
established  from  these  results  why  the  shock  is  not  immediately  plane 
and  discontinuous,  but  it  might  be  due  to  either  an  initial  shock  curva¬ 
ture  because  of  the  bowed  diaphragm  (in  which  case  there  would  be  oblique 
reflections  from  the  walls  trailing  the  shock)  or  due  to  the  finite 
time  required  for  the  shattered  pieces  of  diaphragm  to  spread  out  from 
the  point  of  rupture. 

The  next  assumption  made  was  that  the  temperature  contact  dis¬ 
continuity  represented  no  pressure  (that  is,  velocity)  increment.  The 
record  in  figure  10(d)  shows  that  no  pressure  change  is  recorded  at  the 
point  (about  1.75  milliseconds  after  the  compression)  where  the  dis¬ 
continuity  passes  the  crystal. 

Finally,  the  theoretical  method  is  based  on  the  assumption  that  the 
diaphragm  instantaneously  disappears  such  that  not  only  are  plane  discon¬ 
tinuous  compression  and  expansion  fronts  generated  but  a  plane  temperature 
contact  discontinuity  is  also  generated.  The  schlieren  photographs  and 
pressure— time  records,  however,  show  that  this  assumption  is  not  the  case 
experimentally.  The  pressure— time  records  show  that  the  expansion  zones 
are  slightly  thicker  than  expected  from  theory  (fig.  10),  and  the 
schlieren  photographs  and  pressure— time  records  show  that  the  temperature 
contact  -discontinuities  appear  to  be  a  thick  region  of  mixing  in  the  flow. 
Figure  l4,  in  which  the  reflected  shock  is  being  propagated  through  the 
temperature  contact  discontinuity,  shows  that  some  distortion  of  the 
shock  is  evident,  probably'  a  result  of  the  unequal  temperature  distribution. 
From  figure  10(a),  a  gradual  rise  occurs  where  a  sharp  compression  was 
expected  theoretically  as  a  result  of  the  interaction  of  this  reflected 
shock  with  the.  temperature  contact  discontinuity. 

A  possible  explanation  for  these  results  might  be  that  in  the  case 
of  an  expansion  which  is  not  initially  plane  because  of  the  diaphragm 
bowing,  oblique  reflections  of  the  expansion  from  the  wall  would,  not 
theoretically  catch  up  with  the  expansion.  Also  the  temperature  contact 
discontinuity,  when  not  initially  plane,  would  become  a  vortex  in  the 
flow  because  of  its  entropy  variation. 

The  pressure  records  have  previously  shown  that  the  reflection  of 
the  shock  from  the  end  plates  gave  rise  to  mechanical  oscillation  of  the 
end  plates.  The  pressure  waves  generated  as  a  result  of  this  oscillation 
can  be  seen  in  figure  14  to  be  trailing  the  reflected  shocks. 


CONCLUDING  KEMAKKS 


Experimental  pressure— time  measurements,  shock— speed  measurements, 
and  -schlieren  photographs  have  "been  presented  and  the  results  indicate 
that  the  unsteady  flow  generated  experimentally  conforms  accurately  with 
the  flow  defined  by  one-dimensional,  nonviscous,  unsteady— flow  theory. 

The  agreement  of  experimental  measurements  with  theoretical  calculations 
was  considered  excellent  for  flows  involving  pressure  disturbances 
(including  compression  shocks  or  expansion  zones),  temperature— contact 
discontinuities,  and  reflections  or  interactions  of  these  quantities. 
Slight  discrepancies  were  observed  in  the  correlation  of  the  data  and  were 
believed  to  be  introduced  by  failure  to  generate  exactly  the  idealized 
type  of  flow  assumed  in  the  theory.  The  successful  recording  of  the 
unsteady  flow  properties  was  largely  attributed  to  the  special  equipment 
developed  for  this  purpose. 


Langley  Aeronautical  Laboratory 

National  Advisoiy  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Ya.,  March  17,  19^9 
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APPENDIX 

IN  STR13MENTATI0N 

Apparatus  for  Measuring  and  Recording  Pressures 


Crystal ,  preamplifier,  and  filter The  theoretical  shape  of  the 
pressure  history  to  he  recorded  is  shown  in  figure  10.  The  steepness  of 
the  leading  edge  and  the  flatness  of  the  constant— pressure  portion  of 
the  record  are  dependent  upon  the  response  characteristics  of  the  crystal 
and  preamplifier  combination.  The  crystal  and  preamplifier  combination 
required  to  reproduce  the  pressure  history  shown  is  essentially  one  with 
a  very  low -frequency  response  as  well  as  a  high-frequency  response. 

The  pressure-measuring  element  used  was  a  Massa  M-101  piezoelectric 
pressure  pick-up.  The  first  records  taken  showed  that  the  high-frequency 
response  of  the  crystal  was  sufficient  to  reproduce  satisfactorily  the 
steep  leading  edge  but  that  poor  low— frequency  response  resulted  in  a 
noticeable  drop  in  the  output  during  the  time  of  constant  pressure.  The 
open— grid  cathode— follower  preamplifier  circuit  shown  in  figure  15  was 
developed  to  present  a  maximum  impedance  to  the  crystal  unit  in  order  to 
obtain  low-frequency  response.  The  addition  of  the  100-micromicrofarad 
condenser  across  the  input  further  increased  the  low— frequency  response 
without  noticeably  affecting  the  steepness  of  the  wave  front. 

Figure  1 6  shows  the  input-output  voltage  of  the  preamplifier  alone. 
Negative  input  voltages  are  seen  to  have  resulted  in  distortion  which 
made  it  necessary  to  bias  the  grid  positively  when  recording  calibrating 
voltages  if  a  sine-wave  voltage  was  used.  The  grid- to- cathode  resistance 
is  shown  to  be  effectively  infinite  until  the  grid  is  14  volts  positive 
so  that  leakage  of  the  piezoelectric  charge  was  due  only  to  the  internal 
resistance  of  the  crystal  itself.  Records  of  crystal  leakage  were  made 
by  impressing  a  direct-current  voltage  momentarily  across  the  crystal 
when  in  combination  with  the  preamplifier  and  100-micromicrofarad  con¬ 
denser  and  by  recording  the  exponential  decay  of  the  voltage  output  of 
the  preamplifier;  the  procedure  was  then  repeated  without  the  crystal  in 
the  circuit.  The  time  constant  of  the  crystal  was  found  to  be  55  milli¬ 
seconds;  that  is,  two-thirds  of  the  charge  on  the  crystal  was  dissipated 
in  55  milliseconds,  about  ll-  percent  per  millisecond.  In  order  to  maintain 
this  time  constant,  the  crystal  must  be  dehydrated  after  each  day*s  use 
by  storing  it  in  a  flask  containing  a  dessicant.  The  slope  of  the  trace 
during  the  constant— pressure  portion  of  the  records  shown  in  figure  10 
illustrates  the  very  small  leakage  occurring  during  this  time. 

-  The  appreciable  amplitude  of  hash  observed  on  the  record  shown  in 
figure  17(a)  was  generated  by  the  crystal*  s  resonating  at  its  natural 
frequency  of  about  47  kilocycles  after  being  excited  by  the  steep  wave 
front.  Satisfactory  attenuation  of  the  resonant  frequency  was  obtained 
by  the  use  Of  filters  as  is  shown  by  the  record  of  figure  17(b)  which 
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was  taken  under  conditions  similar  to  those  for  figure  17(a)  except  that 
a  filter  was  used.  The  filter  was  constructed  in  accordance  with  the 
diagram  shown  in  figure  l8  with  parts  that  were  accurate  to  within  ±1  per¬ 
cent  and  functioned  satisfactorily  with  the  different  crystals  used 
during  the  investigation.  With  the  filter  in  the  circuit  some  attenuation 
of  frequency  components  above  5000  cycles  can  he  expected;  however,  the 
rise  time  at  the  leading  edge  would  he  the  only  part  affected  in  these 
records. 

Because  the  preamplifier  output  for  the  range  of  pressures  covered  in 
these  tests  was  of  the  order  of  10  volts,  amplification,  in  addition  to 
that  of  the  oscilloscope  amplifier,  was  not  needed. 

Recording  apparatus.—  At  first,  records  were  made  by  photographing 
a  four-beam  commercially  constructed  oscilloscope.  The  four  beams  made 

it  possible  to  record  simultaneously  the  pressure  history,  a  —  second 

timer,  a  1000  cycle  calibrating  voltage,  and  a  reference  line.  Because 
defocusing  at  high  intensities  masked  many  details  in  the  records  made 
with  the  four-beam  oscilloscope,  it  was  replaced  with  a  single-beam 
oscilloscope  of  higher  resolution  that  gave  records  of  more  detail  and? 
higher  contrast.  When  this  single— beam  oscilloscope  was  used,  the 
reference  line  and  calibrating  voltages  were  recorded  immediately  preceding 
a  test  on  the  same  film  record.  Calibration  indicated  that  the  spot 
deflection  on  the  oscilloscope  varied  linearly  with  input  voltage  to  the 
preamplifier  until  saturation  of  the  preamplifier. 

Moving  film  and  no  sweep  of  the  spot,  rather  than  single  sweep  with 
stationary  film,  was  chosen  as  the  method  of  photographing  the  record. 
Records  were  taken  with  an  NACA  drum  camera  using  an  f/2.0  lens  and 

orthochromatic  film  2=5-— inches  wide.  This  camera  consists  of  a  drum 

lo 

16^  inches  in  circumference  rotated  at  3^00  rpm  ^7  means  of  a  synchronous 

motor  so  that  the  film  speed  is  0.95  1  percent  inch  per  millisecond. 

The  shutter  was  set  for  l/50  second  and  was  synchronized  with  the  drum 
of  the  camera  and  solenoid  trigger  on  the  experimental  apparatus  as 
described  in  the  following  section. 

Synchronizing  circuit.—  Figure  19  shows  a  block  diagram  of  the 
synchronizing  circuit.  A  small  bar  magnet  cemented  on  the  film  drum 
generated  a  trigger  signal  as  it  passed  a  coil  mounted  on  the  drumr-case 
lid.  The  trigger  signal  caused  the  shutter  solenoid  to  be  energized 
which  in  turn  released  the  shutter.  After  a  delay  of  a  few  milliseconds, 
the  external  solenoid  for  actuating  the  diaphragm-rupturing  trigger  on 
the  experimental  apparatus  was  also  energized.  Film  was  loaded  on  the 
drum  so  that  it  began  and  ended  at  the  bar  magnet,  thereby  establishing 
a  definite  time  relation  between  the  trigger  signal  and  the  instant  the 
beginning  of  the  record  passed  the  lens.  The  shutter  was  set  for 
l/50  second  or  one  complete  revolution  of  the  film  drum. 
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The  pulse-delay  circuit  (fig-  20)  which  synchronized  the  schlieren 
apparatus  with  the  diaphragm  "burst  was  actuated  by  a  crystal  near  the 
diaphragm  and  was  set  to  give  the  delay  required  for  the  shock  to  travel 
from  the  crystal  to  the  glass  section  of  the  experimental  apparatus. 

The  delay  was  set  by  means  of  a  10— turn  helical  potentiometer  which  was 
calibrated  by  use  of  a  Potter  interval  timer.  It  was  determined  that 
the  delay  interval  could  be  reset  within  ±10  microseconds  consistently. 
The  pulse-delay  circuit  was  connected  to  the  socket  of  existing  spark 
equipment  used  for  schlieren  photography  so  that  the  output  voltage 
pulse  fired  a  thyratron  which  was  in  series  with  the  primary  winding  of 
an  ignition  coil.  The  ignition  coil  in  turn  fired  a  series  pilot  gap 
which  -set  off  the  main  spark  gap. 


Apparatus  for  Recording  Shock  Speed 

The  speed  of  the  shock  was  determined  by  measuring  the  interval 
between  the  time  the  shock  passed  the  window  of  the  experimental  apparatus 
and  the  time  it  reached  the  end  plate  by  means  of  the  following  arrange¬ 
ment:  A  sine  wave  of  10,000  cycles  per  second  was  impressed  on  the 
vertical  axis  of  an  oscilloscope  to  provide  a  time  scale,  and  two  pulses 
representing  these  two  positions  of  the  shock  were  applied  to  the 
horizontal  scale.  Cfae  of  these  pulses  was  the  variation  in  the  output 
of  a  931  photomultiplier  tube  caused  by  the  bending  by  the  shock  wave  of 
a  beam  of  light  which  emerged  from  a  slit  0.020— inch  wide  and  1— inch 
high,  passed  through  the  windows  of  the  experimental  apparatus  at  a 

point  19^  inches  from  the  end  of  the  low-pressure  chamber,  and  activated 

the  photomultiplier  tube.  (The  light  from  the  slit  was  cut  off  from  the 
photomultiplier  by  a  knife-edge  when  no  shock  was  present.)  The  other 
pulse  was  the  output  of  a  crystal  caused  by  the  impact  of  the  shock  on 

the  end  plate  of  the  72^-  -inch  low-pressure  chamber  to  which  the  crystal 

was  attached.  The  trace  on  the  oscilloscope  was  photographed  with  the 
drum  camera  synchronized  to  the  diaphragm  burst  as  previously  described. 
The  sine-wave  frequency  was  accurate  to  within  one  part  in  100,000  and 
the  outputs  of  the  photocircuit  and  crystal  (without  preamplifier  and 
filter)  were  sharp  enough  for  the  departure  of  the  spot  from  the  sine 
pattern  to  be  read  to  about  1  microsecond  in  1200.  Six  records  were 
taken  at  each  point  shown  in  figure  2  and  the  maximum  scatter  of  the 
records  at  each  point  was  less  than  ±1  microsecond  in  about  1200. 
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Value  of  pg  /p^  calculated  from 
equation  (5)  and  measurements  of 


Figure  2.—  Theoretical  and  experimental  pressure  ratio  across  compression  shock  "before  reflections 
as  function  of  pressure  ratio  across  diaphragm  before  bursting.  &q  =  ai. 
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Denotes  these  quantities  constant 
throughout  an  enclosed  area 


Figure  4.—  Typical  diagram  of  x  against  t  showing  position  histories 
of  unsteady  flow  quantities  following  diaphragm  hurst  in  constant 
area  passages,  slq  =  a^_. 
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Figure  5.—  Pressure  distribution  along  axis  of  propagation  of  expansion  zone. 


Figure  7.—  Sketch  of  apparatus  used  to  study  flow. 
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Figure  8.—  Apparatus  used  to  study  flow. 
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Distance,  in. 


(  c )  72^  —  inch  low— pr © s sure  chauib er . 


Figure  9*~“  Concluded. 
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(c)  —  =  2.0;  72^— inch  lew— pressure  chamber  at  end  plate. 

Pi 


(d)  —  =  2.5;  2k~—  inch  lew— pressure  chamber  at  end  plate. 

Pi  ^ 


Figure  11.—  Typical  pressure  records  obtained  with  standard  test  configu¬ 
ration  except  where  noted.  _  _ _ 

n'aca.  ’ 

L-5981+3 


MCA  TN  I9O3 


47 


(e) 


—  =  1.3;  2^- inch  low-pressure  chamber  at  end  plate. 
P1  d 


(f) 


—  =  2.0;  48-inch  low-pressure  chamber;  crystal  located  in  the 
P1 

side  24  inches  from  diaphragm  in  low-pressure  chamber. 


=  1.3;  48-inch  low-pressure  chamber;  crystal  located  in 
side  24  inches  from  diaphragm  in  low-pressure  chamber. 


the 


=  2.0;  7 2^-  -inch  low-pressure  chamber;  crystal  located  in  the 
side  7|-  inches  from  diaphragm  in  low— pressure  chamber. 


Figure  11.—  Continued. 


l-5981 A 


(i)  —  _  2.0;  7^-— inch  low-pressure  chamber;  crystal  located  in  the 

P1 

side  20  inches  from  diaphragm  in  high-pres3ure  chamber. 


(j)  _  2.0;  72i -  inch  low-pressure  chamber;  crystal  located  at  end 

P1  d 

plate  of  high-pressure  chamber. 


00  52 

pl 


1.3;  72^— inch  low-pressure  chamber;  crystal  located  at  end 


plate  of  high-pressure  chamber. 
Figure  11.—  Concluded. 
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Figure  12.—  Voltage  output  of  crystal  unit  as  a  function  of  pressure 
difference  experienced  at  crystal  as  a  result  of  passing  disturbance. 
Ap  calculated  from  theory.  (Symbols  designating  test  conditions 
are  identified  in  table  II.) 
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(a)  —  =  1.6;  24^—  inch  low-pressure  chamber; 

Pi  2 

first  window  from  diaphragm;  time, 

0.5  millisecond. 


(h)  —  =  2.0;  2k~r  —  inch  low— pressure  chamber; 

Pi  2 

first  window  from  diaphragm;  time, 

0.5  millisecond. 


(c)  ^  =  2.0;  24^  — inch  low-pressure  chamber; 

second  window  from  diaphragm;  time, 

0.97  millisecond. 


Figure  14.—  Schlieren  ^photographs  taken  at  various  times  following  "burst 
of  diaphragm.  Positive  x-direction  is  to  the  Mgtefe-. 
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(d)  =  2.0;  24^— inch  low— pressure  chamber; 

second  window  from  diaphragm;  time, 

1.08  milliseconds. 


p  p.  1 

(e)  —  =  2.0;  inch  low— pressure  chamber; 

P1  2 

second  window  from  diaphragm;  time, 

I.98  milliseconds;  shock  reflected. 


(f)  —  =  2.0;  2^—  inch  low-pressure  chamber; 

Pi  2 

first  window  from  diaphragm;  time, 

2.0  milliseconds. 

Figure  l4.-  Continued. 


(g) 


£o 

Pi 


=  2.0;  24^  —  inch  low- 


pressure  chamber 


first  window  from  diaphragm;  time, 
2.7  milliseconds;  shock  reflected! 


(h)  “  -  2.0;  72“  inch  low-pressure  chamber; 

second  window  from  diaphragm;  time 
4.22  milliseconds. 


(i)  — ' ~  -  2.0;  72^— inch  low-pressure  chamber 

P1  ^ 

second  window  from  diaphragm;  time, 

5.3  milliseconds;  shock  reflected. 


Figure  l4.—  Concluded. 


Figure  15.—  Preamplifier  circuit 
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(a)  Record  taken  without  filter  In  circuit. 


(b)  Record  taken  with  filter  in  circuit. 

Figure  17.—  Experimental  records  showing  attenuation  of  natural  frequency 

of  crystal  by  47-kilocycle  filter. 
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Bar  magnet  cemented 
to  film  drum 
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Figure  19.—  Highr-speed  camera  synchronizing  unit 


